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We define and compute from data the strong couplings of the X(3872) with both of the possible 



quantum numbers assignments J 



1 +,2 . We use these to compute cross sections for J/ip 



resonance scattering into DD*. As an application of the results obtained we revise the calculation 
of the J/4> absorption in a hot hadron gas to confront with recent RHIC observations in Au-Au 
collisions. 

PACS: 12.39.-x, 12.39.Mk, 12.38.Mh 
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I. INTRODUCTION 

Recently the BaBar collaboration took issue with the statement that the X(3872) is a 1"'""'" resonance, as 
was widely accepted, raising the hypothesis of 2 quantum numbers yet to be confirmed by Belle or LHC 
experiments. 

In consideration of this we have discussed in a previous paper ||2] the consequences of the quantum numbers 
assignment suggested by BaBar confronting the expected mass of the 2 standard charmonium with the one 
measured for the X. We also raised the problem of the small X prompt production cross-section predicted 
for a 1^Z?2 charmonium when confronted with the measured one. Indeed, assuming the 2 ^ hypothesis, the 
molecular interpretation of the X(3872) (at odds with data also in the l"*""*" case, as discussed in ^, i^] and 
contested in is ruled out and the charmonium interpretation comes back into play. Also the tetraquark is 
still an open option f^. 

In this paper we will expand on the consequences of the J^*^ assignment studying the decay modes of the 
X(3872) under the hypothesis that it is a 1++ state, let us call it Xi, or a 2 ^ state, call it X2, but making no 
hypotheses on its structure (charmonium, molecule, tetraquark). 

As a first step, we will define a general parameterization of the transition matrix elements, describing the 
known decays in terms of a set of strong coupling constants. Using the experimental information available and 
summarized in Table U as discussed in 6], we will determine the strong couplings in the 1"^"*" and 2 cases by 
the explicit computation of decay widths. In our fits we will also use the data on A" J/ip w decays reported 
inQ. 



X 



Decay Mode 


B 


a{B) 


X -> J/lp TT + TV- 


0.055 


0.020 


X J/iP TT+Tr^Tr" 


0.045 


0.030 


X J/i! 7 


0.0135 


0.0060 


X -> D^D"* 


0.67 


0.13 



TABLE I: Branching ratios B and one sigma errors (j{B) for the observed decays of X(3872) [^. 



We will not attempt any theoretical determination of the strong coupling constants we are going to define. 
This would require to formulate some hypotheses on the structure of the X and the use of approaches such as 
quark models or QCD sum rules. This work could be done elsewhere and confronted with the coupling strengths 
found here. 

We also confirm that using data in [l[ , the negative parity assignment for the X is indeed favored: as opposite 
to an earlier analysis by CDF Q on the J/tj} tt+tt" angular distribution, indicating that both the l"'"^ and 2 
assignments are equally possible, we will show that the 2 ^ assignment would be the preferred one. 

Once we have extracted from data the strong couplings of the Xi^2 to the w and p vector mesons, we will 
calculate the J /ij} {p,uj) — Xi^2 — ^ D^D^* cross-section. This calculation could be of relevance to the study of 
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a classic background to the J /ip suppression signal in heavy ion collisions: the J /ip can be absorbed by a hot 
gas of hadrons including pions and the lighter vector resonances at a similar rate as it is supposed to occur in 
a phase of deconfined quarks and gluons. If one finds that this mechanism is effective at temperatures smaller 
than the Hagedorn temperature, considered as the limiting one for the hadronic matter (see e.g. 0), then the 
J /if) suppression signal should be regarded as the less compelling among the many indications available of a 
new phase of matter produced in heavy ion collisions. 

We will show that only the X2 would contribute effectively to the J/ij} {p,u!) D^D^* absorption, whereas 
the Xi has no significant contribution. It is clear that given the narrowness of the X(3872), the process of 
converting J /ij} {p, ui) into D^D^* due to an intermediate X is too slow to have a mean free path in a gas of p's 
and oj's that is likely larger than the typical size expected for the fireball generated in a heavy ion collision. 

Most of the studies on the hadronic suppression of J/V' are made under the hypothesis that the in-medium D 
and D* mesons have a reduced mass and a larger width [9l-[l7|. This in turn would compensate for the narrow 
width of the X, because of a larger D^D^* phase space - we find that mass variation of the open charm mesons 
allows a sizable contribution of the X2 in the hadronic J suppression. It has been claimed that if one sets a 
definite temperature at which the D and D* masses drop down, one could reproduce the dip observed in Jjip 
production data in correspondence of a critical centrality of the nucleus-nucleus collision [3| • As discussed 
in the section on the comparison to data on J suppression at RHIC, we do not confirm this pattern in our 
analysis. 

We will add the contribution of the X2 to former results on the non-resonant hadronic J / ip suppression due 
to a hot n/rj, p/uj, (j), k'^ ' gas and show its relative weight with respect to that. This study is intended to be 
suggestive of the fact that the XYZ resonances could be of relevance as intermediate states in a number of 
physical processes, as proposed for example in ^2d\. 

In Section [n] we define the strong coupling constants for Xi and X2 and calculate their values using available 
data on the X{3872) decay modes. In Section fill Al using the results obtained in Section |lll we compute the 
cross-sections for processes like J/ip{p,oj) — )■ Xi^2 — >■ D^D°*. Section UlI CI is devoted to the computation of 
the average absorption length of the J/ip in a hadron fireball, taking into account the in-medium properties 
of open charm mesons described in Section IIIIBI We compare these results with those obtained in absence 
of effects modifying the mass and width of D mesons. We include in our analysis resonant and non-resonant 
contributions neglecting their interference. Finally in Section IIIIDI we compare our predictions to the RHIC 
data on Au-Au collisions. 



II. X DECAYS 

We start with the parameterization of the transition matrix elements for the decay processes in Table U in 
terms of coupling strengths whose numerical values are then extracted by comparison with experimental data. 
In the next subsection we discuss the Mfi matrix elements which are related to the transition matrix elements 
Tf^ by 

/ 

As for the normalization of states in Mfi the standard l/v2EV is used. 

A. Transition matrix elements 

We require that strong transition matrix elements are parity even Lorentz scalars obtained by combining the 
momenta and polarizations of the initial and final particles. The conservation of angular momentum fixes the 
decay wave of A — >■ BC: J a = {Jb ® Jc) ® ^bc, ^bc being the relative orbital angular momentum in the final 
state. For each unit of orbital angular momentum in the final state there must be factor of a spatial component 
of the momentum in the transition matrix element. Here and in the following we use the notation ip and J/ip 
interchangeably. 

The J^'-' — 1++ case. The decay X — ?> ipV, with V — p,uj is a.n £ ~ decay, since from the point of view 
of the quantum numbers it corresponds to 1+ — > 1^1^. There is only one combination of momenta and 
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polarizations which has all the properties we enumerated above ^ 

We,p)V{r^,q)\XiX,P))^g,^,. e^'^P^ A,,(P) eUp) r^;{q) P,. (2) 

The decay X — > D^D^* is also an £ = process, since it corresponds to 1+ 0^1^. The matrix element can 
be written in terms of a second coupling strength, g^^u, , as follows 

{D\p)&*{e,q)\X{\P))^g^^^. \^{P) e^iq). (3) 

In order to conserve charge conjugation one should consider the final state D^D^* + D'^D^* . As explained in 
Sec. ID ll of the Appendix, we can consider only the D^D^* component of the final state in what follows. 

The J^*^ — 2 ^ case. In this case, both the decays X — )■ '^}:V and X —J- D^D^* are I — 1 processes, since 
they correspond to 2^ — > 1^1^ and 2^ — >■ 0~1~ transitions respectively. 

The spin of the X is described by a symmetric traceless polarization tensor tt'^" satisfying P^n^^ = 0. In 
the rest frame the five independent components can be set in a 3 x 3 traceless tensor tt'-' . For the sum over 
polarizations we have (2l| 

^7r^^(fc)7r*^(fc) = -^{gt^agiylS + gti/Sgi^a - gfiuga/^) - {9f^aKkp + gi^pk^ka + gf_ipKka + (4) 

pol 

12 2 

+ -^{gfj-u + ^k^k^){gai3 ~i i^kakp), 

with fc^ = m^. 

For the decay X — > -01/, we have to determine the transition matrix element {'il){e^p)V{rj,q)\X{'K^P)). There 
are two ways of combining momenta and polarizations which give a parity even Lorentz scalar^: 

{i) a polarization vector contracts with the left index of the tt tensor. If that of the ip we have 

e*"(p)7r„^(P)e^^''>.9pr;;(g) (5) 

or if the V one does 

r,*^{p)n^^{P)tP''P'^ q^Ppeliq)- (6) 

(m) a momentum contracts with the left index of the vr tensor. One can have 

p" 7:^,{P) e^'"'^ q, e;{p) ^:{q) (7) 
and the remaining combinations of momenta obtained by replacing pq by pp, qq, qp. 

Since we have a P - wave decay, we should not have non-zero terms proportional to Piqj, where i and j are 
spatial indices. In the X rest frame such terms are absent in Eq. ([5]) and ([6|) since they would be proportional 
to p X q - which vanishes only in the rest frame of the decaying particle. The only non zero combination of the 
type of Eq. ([7]), not containing the piqj terms is 

Q"7r„^(P)e^''''^P.e;(p)r/:(g), (8) 

where Q — p ^ q and P = p + q. 

In conclusion we find that there are only three invariant amplitudes one can form by combining these tensors 

Ti = (p) t:^^ (P) e'^'^P" p, q^ < {q) + ,7*" (g) (F) e^'^"- p, q, C (p) , (9) 



^ In the rest frame of the decaying particle Pc = {rri-^ , 0) and one can write 

{^ie,p}Virj,q}\XiX,P)) = g,^^ A,(P) = g,^^, (A(P)xe*(p)) -ri'iq) 

which is the scalar product of two polar vectors, the first coming from the vector product between an axial vector A and a polar 
vector e. Moreover the above expression does not contain any spatial component of the momenta and thus accounts for an 
S - wave process. 

^ An e'^'^P" tensor is needed to obtain even parity. Moreover one cannot contract the two indices of the symmetric tt tensor with 
two of the indices of the completely antisymmetric e tensor. 
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Ts = Q" n^^iP) e'^-P" e;{p) tjUq)- (11) 

which carry three imphcit polarization indices. The first two correspond to the sum and the difference of Eq. ([5]) 
and ([6]), which it turns out to be useful to further reduce the number of independent tensors. Indeed one can 
show that Ti and T3 are one and the same tensor. To do this we prove that the following relation among sums 
over polarizations holds 

(^TiT*)^ = (^|rin(^|T3H. (12) 

pol pol pol 

The above condition implies that the two tensors are equal up to a constant if the sum over polarizations has 
the properties of an inner product. The Schwarz inequality states indeed that for all vectors v, w 

\{v,w)\^ <{v,v){w,w), (13) 

where (•, •) is an inner product: the equality holds only if the two vectors are linearly dependent, i.e., if they 

are parallel. Given two vectors v, w £ C", the inner product is defined as {v,w) = X]^=i ^^n^n = ^n=i "^n^n- 
Here we are evaluating sums over polarizations, labeled by n, which means 

5]T,i;* = 'f]%(")(7;("V, (14) 

pol 



where we are summing over the 5 polarizations of the X and the 3 of the vectors. Therefore Eq. (IT^ implies 
that for each polarization configuration, Ti and T3 are equal up to a constant, and we can choose one out of 
the two for our basis of linearly independent tensors. We choose to keep T3 and eliminate Ti. The final choice 
for the parameterization is 

{^{e,p)V{v,q)\X{n,P)) ^ g,^,T2+ g'^^^Ts, (15) 

where V = p,u!. 

Finally we consider X — )■ D^D^* . One can easily build a parity even Lorentz scalar by contracting the tt 
tensor with the D^* polarization vector and the momentum^ 

= .9,,„.^^''e;(g)p.. (16) 

B. Decay widths and determination of the strong couplings 

The J^^ = 1++ case. Since ui and p have different isospin quantum numbers in principle one needs to use 
different couplings to describe these decays: g^^^ and g^^^. To determine these two values we write the partial 
decay widths for X J/i/j p and X — > J /if) oj as in Eq. (|B10|) and (jClSp in Appendix [B] and [C] respectively. 

For X — > J/-0 p — > J/i/' TT+TT^ we have 

1 1 /" j„ Y^i/„/,„^„MT^\i2„*^™2 ™2 „^l mpFp jg(p -» tttt) rup p*{s,ml+,ml_) 



(17) 



^ The even parity can be easily understood. In the rest frame of the X one has 

(DO(p)DO*(E,<7)|X(7r, P)> = g^^^, 7r'^e*(g)p, = g^^^, [a ■ e*{q)) (X) {v ■ p) 
where a and v are an axial and a polar vector respectively defined by tt'-' = (g) + ® . 
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Here B denotes a branching fraction, Tp is the width of the p resonance, {'ipp{s)\X) is the transition amphtude 
of the previous section ^, and p* is the decay momentum in the X rest frame, given by 



P {x,y,z) = — , (18) 

where the Kallen function is 

A(a;, y, z) — x'^ + y'^ + z'^ — 2xy — 2yz — 2xz. (19) 

In the calculations we will substitute mpTp — > -^s Vpis) — > {s/mp)T p, the comoving width (see Eq. (|B17p in 
Appendix [B|) . Similarly for X — > J/V' lo ^ J tt^tt^tt^ we get 

r(X ^ J/V' tt+tt-tt") = I dsY,\{i^oj{s)\X)\^p*{ml,ml,s) 



pol 



1 m^T^ B{uj ^ 3tt) <i>(^)(^, TO^+,TO^-,m^o) 
TT (s - )2 + (m^r^)2 $(3)(TOj^,m^+,m^-,m^o) ' 



(20) 



where m^Ti^ — ^ (s/m„)rt^. The meaning of is explained in Eq. (jC7p . Appendix C. 

The width of X ^ D°D°* -J> D°D°tt° can be written as the one for X J/^p p J/tp tt+tt". Using the 
expressions for the invariant amplitudes in terms of the couplings constants in the preceding section, we obtain 
g,^^ = 0.14 ± 0.03, 5,,,^ = 0.36 ± 0.01 and g^^^, = (5 ± 1) GeV. 

The J^f^ = 2- "+ case. We will use four different couplings to describe the decays X ^ J/tp p and 
X J/tp Lu: g^^p, g'^^p and .92^^: 92^,^- '^/^ channel, one can assume that the decay proceeds 

through a hadronic channel: X first decays to J/tp p or J/tp uj and later p or uj convert into a photon, using 
vector meson dominance 

{J/tP j\X) = {^\uj)\{J/^ u:{q^ = 0)\X) + (7|p(g2 = 0))^(J/V' p\X) 

ml ml ^^^^ 

= ^{J/tP uj{q' - Q)\X) + ^{J/tP p(/ - 0)\X). 
ml mj, 

We use the decay constants for p and uj derived from the e+e^ partial decay width of the two mesons: /p = 
0.121 GeV^ and ft^ — 0.036 GeV^ (2^. The matrix element for the decay of X — > J/-0 7 is thus also written 
in terms of g^^^, g'^^^ and (/j^^, Sa^^- We are left with four couplings to be determined and only three input 
values for the branching ratios: B{X tpu), B{X tpp) and B{X — tp^). To perform the fit of the coupling 
we therefore use the data on the 37r invariant mass spectrum taken from 1] . 

In [ij 37r events are selected from a sample of J/tp uj events with an invariant mass in the interval 3.8625 GeV < 
mj/^ ^ < 3.8825 GeV. To perform the fit we simulate the decay of a 2 ^ particle extracting its squared 
mass Xi = mf randomly with a Breit-Wigner distribution centered at m-^ = 3.8723 GeV and with a width 
= 0.003 GeV 0. For each value x, we require that x, > Q and that 3.8625 GeV < < 3.8825 GeV. 
Having assigned mf the expected number of 37r events with a definite invariant mass m\^ = s, is proportional 
to the distribution with respect to s of the decay width T{X — > J/tp 7r+7r~7r'') 

, 2 ^ dV{X J/tp tt+tt-ttO) 
iV, {ml^ = s) cx ^ '-g 1 , (22) 

which can be computed using Eq. ()20p . Neglecting the overall numerical normalization we obtain 

N,{ml,^s) ^\Y.\{^J/tl,uj\X{m\))\^- ^ p* {m\,m%,s)'i^^\^s,mX .m' ,ml) (23) 



By (i/'p{s)|X) we mean {'tp{e, p)p{ri, q)\X{X, P)) with = s. s is thus the invariant mass of the tttt pair coming from the p. In 
what follows we will use the same notation for the transition matrix element to a final state containing an unstable particle. 
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if rrii > + ^/s. Thus the total number of events at fixed s is 

N{ml^ ^s)^Yl N^{rnl^ = s) 0{m., - - Vs)- (24) 

i 

In Fig. [1] we show the agreement obtained with data (x^/DOF ~ 4.03/4) and we compare it with the experi- 
mental fit obtained using a Blatt-Weisskopf factor to account for the 1 = 1 decay, as was done in 



14 - 
12 - 
10 - 



0.745 0.750 0.755 0.760 0.765 0.770 0.775 0.780 

m3;.(GeV) 

FIG. 1: Our fit from Eq. (|24|) (Red Triangles) compared with experimental data (Black Disks) and the fit in ^ (Blue 
Diamonds). x^/DOF = 4.03/4 



To compute the normalization factor we exploit the partial decay width of r{X — !• tpuj) = B{X ipuj)rx 
as written in Eq. We obtain g,^^ = (1.58 ± 0.16) GeV"^ and g'^^^ = (-0.74 ± 0.34) GeV^^ Using the 

known experimental data on B{X ipp) and B{X — >■ ipj) we obtain two possible solutions for g^^^^ and g'^^^- 
Since the J/ipp -> X2 — > D^D^* cross-section turns out to be roughly the same using the two sets of couplings, 
we choose to use one of them, namely g,^^ = (-0.29 ± 0.08) GeV"\ g[^^ = (0.28 ± 0.09) GeV"^ 

If one fits the same data set assuming J^'~^ = 1++ a x^/DOF ^ 9/4 is obtained, which means that the 
probability of the 1+"'" hypothesis is smaller by a factor of 6 than the 2 ^ one. For the decay X — > D^D^* we 
use the same method to extract the coupling and we obtain g^Dw ~ ^67 ± 50. 

The results are summarized in Table HIl As a consequence of the fact that B{X — > D'^D'^*) > B{X — > -^w) 
and B{X — > D^D'^*) > B{X ipp) we find that the adimensional coupling g^j^jj, is much larger than g^^^ and 
g-^^p. On the other hand all the dimensional couplings turn out to be of the same order of magnitude of the 
mass scales involved. 



coupling^^^^ 


1++ 


2-+ 


f (.7)Dr>. 


(5 ± 1) GeV 


267 ± 50 




0.36 ±0.01 


1.58 ±0.16 GeV"^ 


-0.74 ±0.34 GeV"^ 




0.14 ±0.03 


(-0.29 ±0.08) GeV"' 


(0.28 ±0.09) GeV"' 



TABLE II: Fitted values for the effective couplings of X(3872) to DD"*, J/ip uj and J/ip p for the two J^'~^ assignments. 



III. AN APPLICATION TO THE J/^|) SUPPRESSION IN HOT HADRONIC MATTER 



Recently the PHENIX collaboration published new data on the J/V' suppression in heavy ion collisions 
observed at RHIC [23j , which have been discussed for example in . These data, together with the upcoming 
ones from the LHC - ALICE collaboration, have encouraged us to consider, as a possible application of the 
determination of the Xi^2 coupling strengths, the study of the contribution of the X(3872) to the J /ip suppression 
by a hot hadron gas and to revise some previous results on this topic. 
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A decrease of the J production in heavy ion coUisions is one of the first quark giuon plasma discovery 
signals suggested in the literature [2^. However processes like 

(7r,?7,p,w,0,x(*\...) ^ (25) 

which are at work in a hypothetical hadron gas formed in place or after the deconfined phase of quarks and 
gluons, may also provide a source of attenuation of J/V' - an antagonist signal to the standard one of quark 
gluon plasma suppression. These contributions might also take place at a different stage of the hadronization 
process - once the plasma has converted into hadrons under the hypothesis that the hadron gas is itself in 
thermal equilibrium. Such situations have been extensively studied in the past. Here we take into account also 
the in-medium effects on the open charm mesons discussed in |9|-[l7| and update analyses such as those in p6l - l32j 
and lattice studies like [ssj. Similar studies can be found in [sjor in where a critical temperature (Mott 

transition) is introduced to reproduce the dip observed in J/tp suppression in correspondence of a particular 
centrality; see also [ssl - fSTj . 

In contrast to previous works, we consider here a resonant channel mediated by X12 which turns out to be 
relevant if in- medium effects on open charm mesons are considered. Due to the narrowness of the X(3872), one 
would expect the contribution of the s - channel processes 

J/ip {p,uj) X{i872) ^ D°D°* (26) 

to be negligible. Nevertheless they can be enhanced because the properties of open charm mesons change when 
propagating inside a hadron medium. In particular their masses are expected to decrease, lowering the Z3°Z)°* 
threshold. Also the non-resonant modes are affected by in-medium D meson properties. Hence we re-analyze 
some results previously obtained [26l - [29l Issj , in particular those in [l^ . 

We will obtain an estimate for the cross-section for the process of Eq. (|26p . using the couplings derived in 
the first part of this paper. In the following we will briefiy review how the properties of open charm mesons are 
expected to be modified inside the hadron medium. Using the results found in the literature we will quantify 
the effect of this dissociation process and update the estimates on the non resonant channels of Eq. (^5)) with 
respect to those given in [l^j ■ We will compare the predictions obtained with the experimental data in the last 
subsection of the paper. 

A. Cross-sections 

The cross-section for the process of Eq. ([26|) . depicted in Fig. [21 reads as follows (see Appendix lA)) 

- 1 1 1 A^/^fs m? ) 

a{J/^p^DD^*) = -— — ^— f[s,ml,mlm\^^,ml,ml,) ' °' °' ^ (27) 

9 2 A^''^(s, m,^, m^) 'v f 1.2 u d iq^j. g 

where /(s) is the integral over the scattering angle of the sum over polarizations of the squared matrix element 
/(s,™V"™p'™^i,2'™D'™D*) = J '^cos6'^ |A^viaXi,2p('S,"i^,"ip,m^^^,m^,m^. ,6*). (28) 

pol 

For the matrix elements we use the couplings reported in Table |TT1 The resulting cross-sections are shown 
in Fig. as functions of Ep or E^:, the energies of the p and the uj in the rest frame of the J/ip: s = 
'^l,uj +^% + '^m^Ep^^. 

The functional behavior of the cross-sections shown can be explained as follows. At small values of the energy 
of the incoming p,uj the "exothermic" peak appears^: the threshold energy of the reaction m^-, -\-m^, is indeed 
smaller than the minimum value of -^/s, namely mp + m^, so that the divergence in the fiux factor is located at 
a larger value than the threshold one. At higher energies, s >> m^, the fiux factor behaves as 1/s , whereas 
the phase space is approximately constant (A^/^(s, 0, 0)/s ~ 1) so that 

a{s) ^ - X f{s) as s >> m\. (29) 

s 

Here comes the difference between the 1++ and the 2 ^ assignments. In the Xi case at high energies /(s) ~ s° 
giving (t(s) ^ 1/.S. If instead X = X2, f(s) ^ giving instead (j{s) ^ s^. 



In Fig. [5] the peak is not resolved because of the x-scale chosen. 
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FIG. 2: Feynman graph for J/V^ (p,w) ^ X(3872) ^ D°D°* 
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FIG. 3: Dissociation cross-section of J/V' into 
open charm mesons mediated by Xi as a func- 
tion of the energy of the p in the rest frame of 
the J/-0 {g,^^ = 0.14, g^^^, = 5 GeV). The low 
energy 'exothermic' peak is present. 



FIG. 4: Dissociation cross-section of J/ij) into 
open charm mesons mediated by Xi as a func- 
tion of the energy of the u in the rest frame of 
the J/^ {g,^^^ = 0.36, g^^^, = 5 GeV). The low 
energy 'exothermic' peak is present. 
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FIG. 5: Dissociation cross-section of J/ijj into 
open charm mesons mediated by X2 as a func- 
tion of the energy of the p in the rest frame of the 
Jl^ (ff2,.p = -0.29 GeV-\ g'^^^ = 0.28 GeV"! 
and ~ 267 ± 50). If one uses the other set 

of couplings for p the cross-section is roughly the 
same. 
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FIG. 6: Dissociation cross-section of J/tp into 
open charm mesons mediated by X2 as a func- 
tion of the energy of the uj in the rest frame of the 
J/^ (52^. = 1-58 GeV-i, g'^^^ = -0.74 GeV"! 
and (?2DD* ~ 267). Consider that we are actually 
concerned only with relatively low energy p and uj 
mesons in a Hagedorn gas. 



8 



B. In- medium properties of open charm mesons 



The modifications of tlie masses and decay widtlis of open cliarm mesons D'^ and D^* inside a hot pion gas 
have been computed for example in [10], following the approach discussed in [sl]. Indeed the presence of a gas 
of light hadrons, such as tt's, can sustain scattering processes which involve _D-mesons leading to a modification 
of their masses and widths. These two quantities are both related to the self-energy diagrams, which can be 
written at finite temperature as the thermal averages of the resonant part of the Dtt^'^ forward scattering 
amplitude. 

The decrease of the mass and the increasing decay width for both the D-mesons found in are shown in 
Fig. [7] and Fig. H 

The authors of [13, [3, obtained similar results but with a different approach. They assume that the shape 
of q'q interaction potentials, responsible for the binding of mesons, is sensitive to the temperature. Thus it can 
happen that some discrete levels, corresponding to different cq excitations, are shifted into the continuous part 
of the spectrum becoming metastable states with different masses and non-vanishing widths. Each Z?-meson 
excitation undergoes this transition at a different critical temperature: AM cx — (T — T^)0{T — T^). 

Since we do not find any relevant differences on the J/ip dissociation by using the two approaches, we will 
consider only the one 



a a ny 
in [13 





FIG. 7: In medium mass modification computed 
in for D° (solid line) and D"* (dashed line). 



FIG. 8: In medium total decay width computed 
in O for D° (solid line) and D°* (dashed line). 
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FIG. 9: Logarithmic derivative of the total decay width of X{3872) as a function of the temperature in the case X = Xi 
and = m^{T) + m^, (T) (dashed line) and X = X2 and = 3872 MeV (solid line). 



The broadening and shifting of the masses of the two open charm mesons lead to a modification of the decay 
width and mass of the X(3872). Since we do not have clues on how in- medium effects would modify the mass 
of a X12 tetraquark, we simply assume that if X{3872) — Xi it is a D^D°* molecule; if X(3872) = X2 it 
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is a charmonium state. We remind here that the 1++ assignment is severely at odds with a 2 "^Pi standard 
charmonium interpretation essentiahy because of the small radiative transition rate X — J/'0 7 with respect 
to what expected. In the molecular interpretation the mass of the Xi is directly related to the sum of the 
masses of the and D^* and thus it will decrease with the temperature. In the charmonium assignment (and 
likely also for tetraquarks) one might expect the mass of the X2 to be almost stable with temperature. This is 
because X2 would be the 1^D2 charmonium radial ground state and Debye screening is not expected to alter 
the lowest lying levels [25;ij. The D^D^* width can be computed as 



X w(si,™„(r),r^(T)) w(s2,m^.(r),r^.(r)), 

where by BW we mean the standard normahzed Breit-Wigner function 

BW{s, m, r) = — — ^ (31) 

and s™'" and s™'*'^ are fixed by the kinematics. We show the results in terms of the logarithmic derivative of 
the total width of X\,2 with respect to the temperature, see Fig. [SI in-medium effects make the X2 become 
much broader than X\ . This fact can be understood by taking into account the dependence of the decay width 
on the masses of the particles in the final state. The phase space volume is enlarged proportionally to the 
decay momentum p* . As for the matrix element, if J^'^ = l"*""*" the X — D^D^* decay has £ = Q and thus 
\M\^ ^ constant, while if X is a 2"+ state it has ^ = 1 so that \M\^ cx p*^. Thus if J^<^ = 1++ then Tx oc p* , 
instead if J^^ = 2"+ then Tx (xp*^. 

To summarize, the fact that the charmonium X2 mass is not affected by the medium makes the X2 D'^D'^* 
P - wave decay much larger because the Z?" and masses are instead sensitively decreased in the finite 
temperature medium. 



C. Comparison to data on J/tl^ suppression at RHIC 

The average absorption length (mean free path) of the J/iJ; due to the presence of a p meson gas at temperature 
T is the inverse of the thermal average of the product of the density number p of p mesons and the cross-section 
a (given in Eq. (pS)) ): 



)t = (2sp-f 1) 



d^Pp a{Ep) ^ 2sp4 
(27r)3 e-f^p/'^^^ - 1 27r2 



dE, 



PpEp(T{Ep) 



(32) 



Here the kinematics imposes that 



£;;^'" = max 



V ■ 



2ml 



(33) 



By numerical inspection we have found that it is safe to cut-off the integrals at E"^^ = 1.5 GeV and E'^^ = 
3.5 GeV for J^*-^ = 1++ and J^'~' = 2 ^ respectively. The difference between the two values for E'^^^'^ can 
be understood noting that the cross-section diminishes as the energy grows if = 1++, while it grows with 
energy for J^'-' =2 ^ . 

Given that the masses of the and D'^* mesons are supposed to change with the temperature we need to 
take into account this effect in the calculation of the thermal averages. We average the absorption length over 
the Breit-Wigner distributions of the D and of the D*: the formula for {p<j)t is therefore 



2s„ 



27r2 



dsy 



dsi 



dE, 



Pp Ep (T{Ep, si,S2, (T), (T)) 



(34) 



X BW{si,m, (T), (T)) BWis2, m^, (T), T^. (T)). 
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As already mentioned in the previous subsection, we report only the results obtained using the masses and 
widths of the D- mesons computed in jTo| . If one uses the discontinuous functions for to^(T) and T^{T) 
proposed in [H, [l^, the values obtained for {pcr)^ are of the same magnitude. Moreover, regardless of whether 
the non-resonant channel is included, {pcr)j, does not show any discontinuity that can help in fitting the observed 
dip in the experimental data, contrarily to what shown in [isl. [l9|. The same holds for the lu. 

In Fig.[Tn]and[Tl]we show the results for the inverse average absorption length for the resonant J/^ suppression 
mediated by X2 and initiated by p and w respectively. For the Xi case we find the effect is negligible, since the 
in-medium Xi is still too narrow for p and to to effectively dissociate the J/ip into open charm mesons. 

We also update the estimates for the non-resonant channels enumerated in Eq. ((25l) as discussed in In 
Sec. ID 21 of the Appendix, we give some details on the counting rules for all the J/ip absorption processes we 
consider in the hadron gas. 

In Table IIIII we give a summary of the results for the inverse mean free paths. The contribution of the Xi 
is negligible whereas the contribution from the X2 resonant channel is of the same size of the sum of the non- 
resonant channels up to about T = 120 MeV. With the growing temperature the resonant contribution is found 
to weight less than the non-resonant ones, reducing to a 15% of the non-resonant total at about the Hagedorn 
temperature T 170 MeV. We remind the reader that we have neglected possible interference between resonant 
and non resonant channels. 




0.10 0.12 0.14 0.16 0.18 0.20 0.10 0.12 0.14 0.16 0.18 0.20 



T(GeV) T(GeV) 



FIG. 10: J/ip p D°D°*: Inverse average ah- FIG. 11: J/ip uj D^D^': Inverse average absorption 

sorption length for X2 hypothesis, using = length for X2 hypothesis, using = 3.8723 GeV. 

3.8723 GeV. 



J/V' 

We now take into account the recent RHIC data on the so called nuclear modification factor i?^^^ , reported 

in (23I I as a function of the number of participants in the collision. The quantity rJ^^ measures the ratio of 
the J/^p yield in A -I- A and pp collisions scaled by the number of nucleon-nucleon collisions. We will consider 
only Au-Au collisions at RHIC, due to their higher statistical significance. We will also reconsider the old data 
on Pb-Pb collisions from NA50 [s^ to show how the picture has changed in the last years. 

Refs. [10, im have also considered the possibility that a recombination mechanism could compensate the 
J/ip suppression due to QGP, making the drop in RHIC data less evident with respect to NA50, where this 
mechanism is expected to be weaker due to the much smaller energies involved. However in [24| it was shown 
that the recombination effects are of the same order of magnitude as the experimental uncertainties and thus 
they can be safely neglected. 

The geometry of the heavy ion collision is shown schematically in Fig. I12[ which depicts the time-evolution 
in the center of mass frame. The impact parameter, b, is defined as the transverse distance of the centers of 
the two nuclei. We consider the J/ip to be created with Feynmans a; ~ 0, during the overlap of the two nuclei. 
These particles have to overcome absorption from the column density of nucleons of extension L. In the center 
of mass frame the length of the column is L/j. In the same frame, the density of nucleons is /0nuci7, so that the 
absorption factor is Lorentz invariant and given by exp (— /OnuciCnucii): see |42| . The nuclear absorption cross 



Cnuci, has been determined in 43| from the behavior of the cross-section for p + A ^ J/ijj + Anything 
and d -|- Au — >■ J/tp + Anything 



nucl 



(3.5 ± 0.2) mb. (35) 
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T (GeV) 


(fm-^) 






(pcr)f' (fm^M 


(pcr)t (fm^M 




(pcr)5J^ (fm^M 


0.150 




0.00700 


0.00182 


0.00244 


0.00052 


0.00469 


0.01648 




0.00435 


0.00801 


0.00212 


0.00268 


0.00052 


0.00580 


0.01908 


0.155 




0.00948 


0.00239 


0.00341 


0.00074 


0.00607 


0.02208 




0.00520 


0.01101 


0.00280 


0.00375 


0.00074 


0.00753 


0.02565 


0.160 




0.01267 


0.00311 


0.00467 


0.00102 


0.00774 


0.02920 




0.00622 


0.01478 


0.00365 


0.00516 


0.00102 


0.00967 


0.03402 


0.165 




0.01672 


0.00398 


0.00631 


0.00138 


0.00977 


0.03817 




0.00737 


0.01959 


0.00470 


0.00670 


0.00138 


0.01224 


0.04456 


0.170 




0.02183 


0.00505 


0.00842 


0.00186 


0.01219 


0.04935 




0.00868 


0.02566 


0.00597 


0.00934 


0.00186 


0.01533 


0.05769 


0.175 




0.02821 


0.00633 


0.01109 


0.00247 


0.01506 


0.06316 




0.01010 


0.03326 


0.00751 


0.01234 


0.00247 


0.01904 


0.07398 


0.180 




0.03610 


0.00786 


0.01445 


0.00324 


0.01845 


0.08010 




0.01175 


0.04270 


0.00935 


0.01612 


0.00324 


0.02341 


0.09400 



TABLE III: Inverse absorption lengths as defined in Eq. (|34p for all the particles in the gas. For each temperature we 
show the results obtained for fixed D-mesons masses (upper entry of each cell) and for decreasing D-mesons masses as 
computed in (lower entry of each cell). Since the decays only into DsDs and we assume that Ds mesons do not 
change their masses and widths inside a hadron medium, the upper and lower entry of each cell are equal. As for the 
resonant contribution due to X2 (first column) we do not report the results with fixed D-mesons masses, since they are 
negligible with respect to the non-resonant ones (NR). We do not consider temperatures higher than the value we use 
for the Hagedorn temperature Th ~ 177 MeV. 

As for NA50 one learns from [3] that 

NA60 , , , , 

=(4.3±0.2)mb. (36) 

For the density of ordinary nuclear matter we take /Onuci = 0.17 fm~'^ [4^. 



(a) (b) (c) 




time 



FIG. 12: Geometry of the collision between two identical heavy nuclei with impact parameter 6. After the two nuclei 
have traversed each other, a thermalized gas of lighter resonances is formed. 

In Fig. HW c) we show the hadron fireball produced by the central collisions of the interacting nucleons [1^ 
(the comoving particles tt, p, lu, ...). The fireball has a transverse dimension, /, approximately equal to the length 
of the overlapping region 

l^2R-b, (37) 

where R is the nuclear radius. The attenuation due to the interactions with the hadrons in the fireball is related 
to the average length that a J/ip has to traverse before leaving it. The RHIC data in ;25j are taken in two 
different rapidity regions: a forward rapidity region 1.2 < \y\ < 2.2 and a mid rapidity region \y\ < 0.35. We 
take, for simplicity, a spherical fireball and we simulate the production of a particle at some point inside the 
sphere and with a given direction of the velocity. Assuming a uniform linear motion inside the fireball, the 
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distance d given the starting point r and the direction of the velocity v, can be written in implicit form as 



\r + dv\ = ^. (38) 

The point on the spherical surface where the particle emerges from the fireball is thus r' — r + dv, from which 
one can compute the rapidity of the J/t/j observed, y rj = — ln(tan(0/2)), with 9 the polar angle associated to 
r' . To obtain the average distance one needs to integrate over the two angles which identify the direction of i) 
and over r, taking into account the constraint on the polar angle of the emersion point given by the experimental 
bounds on the rapidity (yo < \y\ < Vi implies 6q < 9 < 9i). To make our simple simulation more realistic we 
take into account that the distribution of the directions of the velocity is not uniform, but can be approximated 
by f{Vrj.) oc where is the transverse momentum with respect to the beam axis. Finally one has 

SdrSdvf{v^)d{r,v)T{e) 

JdrJdvf{v,)Ti9) ' ^""^^ 

with 

m^l ' ^o<o<o.^ ^^^^ 

I elsewhere 

The result of this computation is d^^^ — 0.4 1 and d^^.^ = 0.3 1 in the forward and mid rapidity region respectively. 
Thus, the attenuation factor due to absorption by the comoving particles is 

X"p,r,?. « exp [-I]i(piCr,)T(;)df„d(.„.d)] : (41) 

the subscript i labels the species of hadrons making up the fireball, pi the number density of the effective (i.e. 
above threshold) particles and (Ti the corresponding J/ip dissociation cross-section. Brackets indicate an average 
over the energy distribution in the fireball. This thermal average is computed at a certain temperature T{1), 
which is given by the centrality of the collision, as we shall explain in detail in the next subsection. 

The NA50 measurements on Pb-Pb collisions were inclusive. Hence one needs to integrate the distance d{r, v) 
over the whole range for the polar angle, obtaining d = 3/8 I, as was done in the previous analysis contained 

As noted before, we can compute the nuclear absorption length, L, as a function of b using NA50 data |47l.l48j 
for Pb-Pb collisions. We report this function in Fig.[T21 Exploiting Eq. (I57)) . one can obtain L as a function of 
/. We can reasonably suppose that the same function L{b) can be used in the analysis of Au-Au collisions at 




FIG. 13: Mean length of the path that sl J/tp pro- 
duced during a Pb-Pb collisions at NA50 must 
travel in nuclear matter as a function of the im- 
pact parameter b [13, |3] . 



FIG. 14: Average number of participant nucleons 
in a Au-Au collision as a function of the impact 
parameter b computed using Eq. (1431) . 



RHIC, since An and Pb have approximately the same radius {Rp^^ =7.1 fm and R^^^ = 7.0 fm). 

Putting it all together, we write the attenuation of the J/ip, due to both comovers and nuclear effects, as a 
function of I according to 
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A 



\l) ^ Co + C X exp[-pnuclO-nucli(0] X <2Xp [-E.^ (pjCri)T(/) ] , 



(42) 



where C is an appropriate normalization constant and Co is an offset. To fit NA50 data we substitute df^^j^ij, 
with 3/81. 

To obtain the experimental data |23l | as a function of I we derive the number of nucleons participating in a 
coUision with impact parameter b from geometrical considerations 



iVpart = 4 



2yl 



(43) 



where A is the atomic mass number and R the nucleus radius, see Fig. [T31 Very similar results are obtained 
using a Glauber Model to relate the number of participants to the impact parameter [39( | . 



D. Hagedorn Gas 

Here we wish to determine the function T{1) to be used in Eq. (pij) . We will describe the fireball as a Hagedorn 
gas of resonances. The partition function of a Hagedorn gas in the classical Boltzmann limit {E >> kbT) can 
be written as [i^ 

ln(Z^') = (^^^ ^ J dm p(TO)m3/2e~'"/^. (44) 

p{m) is the mass spectrum of hadronic states, which has the empirical shape 

= — -^-^e™/^., (45) 

with mo = 0.96 GeV, c = 2.12 GeV^ and = 177 MeV Th is known as the Hagedorn temperature. For a 
recent determination see 0|. As soon as T > the integral in Eq. (|44p diverges, hence this thermodynamical 
description is valid up to T < . Above the Hagedorn temperature the system undergoes a phase transition, 
which can be interpreted as the transition from hadronic matter to QGP. 

From the partition function of Ea. (H^ one can easily obtain the the energy density e(T) 

On the other hand the energy density released in a collision is proportional to the factor 

^^l^mi (47) 

where 

2 arccos (&/2i?) - ib/2R)^l - b^/AK^ ' 
Therefore a simple estimate of the ratio of the energy density for two different values of b is given by 

e(6o) 5(&o/i?) 9{bo/R) ' ^ ^ 

Using the Bjorken relation j46| one can estimate the energy density released in a collision with impact parameter 
b by measuring the transverse energy per unit rapidity 
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where tq, the formation time, is usuahy taken as 1 fm, and 7rr^ is the effective area of the coUision. In ^0] the 
PHENIX collaboration finds that in Au-Au collisions with 90 participants, which corresponds to 6o — 9 fni (see 
Eq. (021)), the energy density amounts to 2.4 GeV/fm'^, thus 

e""""" {bo = 9 fm) = 2.4 GeV/fm^ (51) 
As for the NA50 data on Pb-Pb collisions, we take from [5l| 

e''*™(&o = 9.2 fm) = 1.9 GeV/fm^ (52) 

Using these values we obtain the energy density as a function of the impact parameter e(6). On the other hand 
we know the relation between energy density and temperature e{T) from Eq. ([46]) and thus we can deduce T{b) 
and in turn T{1). We show T{1) for Pb-Pb collisions at NA50 and Au-Au colhsions at RHIC in Fig. [THl It is 
evident that over a wide range of I the temperature is almost constant and below the Hagedorn temperature. 
Now we have all the ingredients to perform a best fit of the experimental data using the attenuation function 
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FIG. 15: Temperature of the Hagedorn gas formed after a Pb-Pb collision at NA50 (left panel) or a Au-Au collision at 
RHIC (right panel) as a function of I — 2R — b, b being the impact parameter. 
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defined in Eq. We show the agreement between experimental data and theoretical prediction in Fig. [12] 

for NA50 data, and in Fig. [JT] and Fig. HH] for RHIC data. 




FIG. 16: Attenuation function for the J/^ yield in Pb-Pb collisions as measured at NA50 (Squares) and as predicted by 
the hadron gas description (Red Line). The best fit is obtained for Co = —0.2 and C — 1.4 giving a x^/DOF = 4.9/9. In 
absence of the resonant contribution from X2 and of the in-medium effects on the D mesons we obtain xVDOF = 5.1/9. 
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2 4 6 8 10 12 

l(fm) 

FIG. 17: Attenuation function for the J/tj) yield 
in Au-Au collisions in the forward rapidity region 
1.2 < \y\ < 2.2 as measured at RHIC (Disks) and 
as predicted by the hadron gas description (Red 
Line). The best fit is obtained for Co = —0.6 
and C = 1.8 giving a xVDOF = 5.5/4. In ab- 
sence of the resonant contribution from X2 and of 
the in-medium effects on the D mesons we obtain 
X'^/DOF = 6/4. In this rapidity region the J/ip is 
reconstructed in the a^^M" mode. 
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FIG. 18: Attenuation function for the J/ip yield 
in Au-Au collisions in the mid rapidity region 
\y\ < 0.35 as measured at RHIC (Disks) and 
as predicted by the hadron gas description (Red 
Line). The best fit is obtained for Co = —0.4 
and C = 1.7 giving a xV^OF = 6/6. In ab- 
sence of the resonant contribution from X2 and of 
the in-medium effects on the D mesons we obtain 
X^/DOF = 6.5/6. In this rapidity region the J/ip 
is reconstructed in the e'^e~ mode. 



IV. DISCUSSION AND CONCLUSIONS 



We have determined from available data the strong coupling constants of the X{3872) under the hypothesis 
that it is a l^~^{Xi) or 2 ^(^2) resonance. The results we find may be confronted with theoretical calculations 
making assumptions on the structure of the X: molecule, tetraquark, standard charmonium. 

We use the coupling strengths we find to explore the potential role of Xi^2 in J/ip absorption processes 
like J/tfj {p,uj) — > Xi^2 D^rP* . Such processes might occur in a hot resonance gas produced in heavy-ion 
collisions. Regardless of the detail mechanism by which the resonance gas is formed, processes as the ones 
mentioned above are mimicking the in-plasma J suppression hypothetically due to the Debye screening of 
the cc confining potential, as first discussed in [2^. Therefore this is a background to the Debye J suppression 
signal. How far can we go with a hadron gas picture in fitting data on J/ip suppression at RHIC? A limitation 
to the hadron gas description might come from the excessively high temperature needed for the gas to account 
for the observed J/ip suppression effect. This was discussed in [29,123]: a hadron gas description fails above the 
critical Hagedorn temperature, the highest temperature for hadron matter. 

The analysis in [2^, [23| was based on NA50 data on Pb-Pb collisions, see Fig. [1^1 where at a centrality of 
about 4 fm in units of I = 2R — b, b being the impact parameter, a drop was observed (actually a one sigma 
effect) in the J/'ip yield in going fr om the three leftmost points towards higher centralities. In the low I region 

< 4 fm) the authors of [2y, [23] used also data on S-U collisions and the approach was to perform a best 
fit in that region (where the hadron gas picture is more reliable, the energy density being smaller) with an 
exponential attenuation function at some temperature T. An unreasonably large T was then needed to fit data 
also at I > 4 fm. Using an Hagedorn gas model the fit was simply very poor. Moreover Ref. ^2] pointed out 
a correlation in the / dependence of the J/ip suppression and the enhancement of strange particle production 
observed in NA50 data. 

Here we take a different approach. We note that the drop at Z = 4 fm observed by NA50 (Fig. [TO]) , is much 
less evident in the recent Au-Au RHIC data in Fig.[T7]and absent in Fig.[THl Therefore we fit the whole data set 
(and not only the I < 4 fm region) with an attenuation function computed in a Hagedorn gas having a limiting 
temperature Th = 177 MeV. As stated in Section IlII CI we are neglecting possible charm recombination effects. 

Actually we find a very good fit to data just using the attenuation functions computed in p7| . This is so 
because the nuclear part of the attenuation function in Au-Au collisions at RHIC is expected to be almost the 
same as that in Pb-Pb collisions at NA50, because Pb and Au nuclei are very similar in size and the J/ip nuclear 
absorption cross sections turns out to be very similar at RHIC and NA50. Moreover to define the dependence 
of the temperature on centrality, we use the energy density produced in the RHIC collisions according to the 
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Bjorken description, and we found it to be almost equal to the one computed for NA50. 

Including some hypothetical in-medium effects on mass decreasing and broadening of open charm mesons 0- 
flTj and a resonant contribution from the J/il) {p,u}) — )■ X2 — >• D'^D'^* channel (about the 15% of the non-resonant 
one) we get slightly larger inverse absorption lengths as shown in Table Hill which altogether slightly improve the 
fit to data, decreasing the x^/DOF from 6/4 to 5.5/4; see Fig. [17] and [TBI In the calculation of J/ip absorption 
we assume that X = X2 has more likely a 2 ^ charmonium interpretation, whereas X ~ Xi has a 1++ molecule 
assignment if only because we have no clues on how a finite temperature hadron medium would alter mass and 
width of a tetraquark particle. In this respect a charmonium X2 gets a much larger width because its mass is 
not modified by the medium, while the masses of its decay products D° and D'^* are. The larger width of X2 
is in turn responsible for the most effective J/ip (p, cj) — >■ X2 conversion which has a very low rate for a narrow 
Xi. The Xi is expected to stay narrow even in medium because: i) its mass lowers as the sum of Z?" and Z?°* 
masses, ii) Xi D'^D'^* is an 5* - wave decay. 

The conclusion which can be drawn from this analysis is that, given the hypotheses we use, data on J/^ 
yield in heavy-ion collisions are likely the less compelling ones in the search of a deconfined quark-gluon state 
of matter because they are affected by a large hadronic background in the sense explained above. 

In order to have a clearer picture, it would be very important to have RHIC data on J/ip suppression for 
a larger number of impact parameter bins, in particular in the intermediate centrality range. It is moreover 
our aim to underscore that X, Y, Z particles being discovered in the last few years might have impact on a 
wide class of elementary processes: we find here that if the X{3872) were confirmed to be a 2 ^ state, under 
certain hypotheses on the behavior of open charm mesons in a hot hadron gas, it would give a non-negligible 
contribution to the hadron J/ip dissociation mechanism. 

Appendix A: J/tJj {p,uj) X(3872) D°D°* cross-section 

We give some details on the formulae used in the text. The differential cross-section for J/ip p D^D^* is 

da{J/iP p ^ D^D^*) = ip^ +P^. - P. -Pp)E |M./,p^.oz,o.r (2^^)X°.„ (2.)^2:;,. ' ^^'^ 

with the flux </> defined by 

(p = ^J{P4, -PpY -m^m2 iy^A(s, ^ m^). (A2) 

We use also 



(27r)32w^ (27r)32a;^. ~ IQtt 
The above formulae leads to Eq. ([?7)). Similarly for w 



Appendix B: X — )^ J/ip p 

Here we report the formulae used for the computation of the width of X — > J/ip p. 

dT{X ^ J/iP TT+vr-) = ^ ^ (Bl) 
2sx + 1 ^m^ 



pol 

where 

^ ' ^ ^ ^ ' {2ttY2E^ {2T:f2Ei {2ttY2E2 

Using the narrow width approximation for the p and the unstable particle propagator 



(B2) 



pol pol "^P' + ^"^"^ P> pol 
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where X^poi K""^"" = dfnn with gp^ a constant number. The phase space factor can be rewritten as 

at". = (2.,. / - - - P. - (5^% W 

The notation d^^^^ is to indicate that we have an intermediate p. Now we observe that 

.0 ^4^(4)/p rf'py, ^(2). /-^ 1 P*(w3„m|,s) 



and 

fr, ^4x(4)f \ ^^Pl ff,(2)^ r ^ 1 p* (a, TO^+ , TO^_ ) 

Thus it resuhs that 



(B6) 



* - 4^ 4^ 7i ■ 

The full decay width is then 

r(X^J/^.+ .-) = ^ ^± / d,yKJ/^,(,)|X)pi^^>^il!!M^ lP*{s,n.l, ml_) 

(B8) 

We can relate g^^ to r(/? — > tttt) by 

9% = 6m^r(p ^ tttt)-— ^ 2 2T 



(s - m2)2 + (mprp)2 ^ p*(m2, , m^) 



pol 

(BIO) 

In the limit of narrow width for the p 



mV 



limr^o?: 2\2~7ZT^ = '^^('5 ~ "^^)- (^^1) 



Eq. (IBlOp is equal to the one we can obtain taking the average of X^poi "^^^ ^) O'^^^ the 

Breit-Wigner distribution of the p meson 

In actual calculations we use mpTp — > ;^rp- Indeed the propagator of an unstable particle A with 4- momentum 
p and which decays to two particles B and C can be written in the form {p^ = s) (ssj : 

1 



where r{p^) = r{A{p^) BC) is 



(B13) 



nA{p-) BC) = (B14) 
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Even if an unstable state cannot be properly put on the mass-shell, the mass of a narrow resonance is still well 
defined and for equal to its mass its decay width is 



T{A{m') ^ BC) = ^'^"f,'"|'"'^V ^("^^'"^B."^!:)- (B15) 



9 9 

167rm 

It is then straightforward to see that 



2 



T^^^2^_ ml g^{p^,m%,ml) \/A(p^ m|, m^^ ^^^^ ^ 



The coupling constant g has the dimension of a mass, to give the right dimension to the width {[g] — M). Now 
let us consider two limits which are relevant to our analysis. The first limit is the one in which both the particles 
in the final state are massless, i.e., much lighter than A, ms ^ mc = 0. In this case the only mass scale of 
the problem is or m^, and thus the only possibility is that g^{p^ , 0, 0) — ap^ and thus g'^{m\, 0, 0) = am\, 
where a is some adimensional constant. The relation (|B16|) reduces to 

r(p2) ^ v^r(m^). (B17) 

ruA 

This relation can be used when the unstable propagating particle has a mass much larger than that of its decay 
products, as in the case of the p. For an interm ediate D* we need to consider the case in which only one of the 
produced particles is massless mc — 0. Since ^A(p2, m^, 0) = {p^ — m\) one obtains: 

r(P^)-^7^4^r(<). (B18) 



{m\ - ml) 



Appendix C: X J/^ u 



Here we report the basic formulae for the computation of the width of X — > J/^ uj 

] 

2^ 



dViX ^ J/^ tt+tt-tt") = \ ^ ^+7r-^°|X)|2d$W, (CI) 

2sx + 1 2m„ ^ 

pol 



where is 

4 

Using narrow width approximation 

pol pol y / 

and we further assume that X^poi K""^^ """^Ip)!^ ^ 5w7r' with gi^^r a constant number. The phase space factor 
can be rewritten as 

<Pp^ TT '^^Pi 



(2^)32^;,' 

where the notation (il>*^^^ is to indicate that we have an intermediate uj. Now 

^ ^ P-\2.)32VIT^(2^)^2i?,-4. m^ ^^^^ 
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and 



i i 

The expression for the three body phase space is the fohowing 

$•^^^(^5,7^1,7712,7713) = y"dw (x+(Vs,7ni, 7772,7773, Cj) - a;_(Vs, 7771, 7n,2,77l3,i:j)) , (C7) 



where 



^""^ "w?"^ "'^ ± ^J{^^ - ml) {oJm i^/s, mi, 7712, 7713) - a;) (^ ^""^^ - cj + cj^ 7771, m2, rris)^ 
a;±(\/s, 7/71, 7)^2, 7773, oj) = 



2v^ + t»Jm (V-5,7ni,7n2, 7773 

with 



(C8) 



/ ^ A "^1 - (W2 + 7773)^ + S 
UJmW mi, 1712, 1713 ) = ■ (C9j 

2vs 



Finally 

1 P*("7i, 777?,,s) 



-'(777x, 777w,, 777jr+, 777^-, 777^0) = — / ds- $^ (VS, 777^+ , 777^- , 777^0 j . ((--lO) 

/TT J AtT 777^ 

The full decay width is then 



p„l -■ --^x (s- 7772)2 + (m^r„)2 1^^^^^ 

X 'l><^^(^/s,7n^+,7n^-,77^^o). 
We can relate to r(tj ^ Stt) 

= 6777„r(w ^ TrTrTr)-^:— -, (C12) 

^y-^) (777^ , 777^+ , 7/7^- , 7/7^0 ) 

thus giving 



pol 

$'^'(yi,7n^+,7n^-,7R^o) 



(C14) 



In the limit of narrow width for the w, Eq. (jC13p is equal to the one we can obtain taking the average of 
EpoiK'^/V' '^(•s)!-^) 7(7-^, s) over the Breit-Wigner distribution of the uj meson (using the comoving 
width of the w) 



r(X^ J/7A^+vr-7rO) = -^-l^ / ^|(J/V.c.(.)|X)|V( 



777v, 7773,,S)- 



1 rn^^i^ B{uj -> tttt) 



(5- 777^)2 + (^^r.)^- 

(C15) 

In actual calculations we use rrii^Ti^ — > -^F^. 



Appendix D: Multiplicity rules 

1. X(3872) D°D°' 

Let us consider the decay X{3872) — > D^D^* . Since X has even charge conjugation, whatever its spin is, the 
final state into open charm mesons needs to be 

^ \D^D^*) + \D^D^*) ^ ^^^^ 
\/2 
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The matrix element is 

(/IX) . 

Assuming that 

{D^D°*\X) = {D°D°*\X), (D3) 
the sum over polarizations of the squared matrix element is 

^|(/|X)p = 2^|(i?°5°*|X)p. (D4) 

pol pol 

When we compute the cross-section for J/-0(p, w) Xj D^D^* we actually consider the transition 
J/tplfyUj) — >■ /. The flavor wave function for the p meson is 

Since the neutral D mesons contain only the u quark {\D'^) = \cu) e \D'^) ~ |cu)) only the uu component will 
contribute to the transition matrix element 

M - ^(/I^P). (D6) 
Summing over polarization the squared matrix element one obtains 

pol pol pol pol 

Thus one can reabsorb the factor 2 into the gj^^, coupling. 

2. Non resonant processes 

We consider the t - channel processes of Eq. (^5]) . We computed the average absorption lengths for each of 
the particles in the initial state, A = tt, 77, p, 0, K(*) using the couplings defined in |27| 

(pcr,//^ a^ddIt = (2sA + 1) y (^2^)3 e^JJ^^T _Y 

Depending on the flavor content of each meson in the initial state one can define the possible open-charm mesons 
configuration in the final state. The flavor wave functions of the mesons we considered are the following (53 | 
(we neglect the ss component of the r] meson, since the contribution of the associated flnal state, Df^*^ D7^*\ 
is small compared to the one coming from the {uu + dd)/^/2 component, i.e., Z30(*)i)°^*' or £)=•=(*) 

7r+(p+) = uJ, T:-{p-) = ud, 7rO(pO) : " 



V2 
uu + dd 



V 



w ~ — — 

(j) ss 

K° = sd, = sd, K+ = us, K~ = us. 

The multiplicity coefficients associated to the possible final states fi for each initial particle A are summarized 
in Table HVl Given these coefficients the total dissociation cross-section for the initial particle A can be written 
as 

^A^Y.''^ ''AJ/^-.U- (DIO) 
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We can summarize all the contributions as follows 

C7^ = 3 X [a(J/i/' TT DD) + 2a{Jltp tt DD*) + a{J/^p n -S- D*D*) , 



<7p+co = 4 X a{J/i; p DD) + 2ct( J/V' p DD*) + a{J/'ip p D*D*) 
ar, = aiJ/ip r] DD) + 2a{J/il> rj DD*) + cr( J/V' V ^ D*D*), 
= <7(J/V </> ^ D;Dt) + 2a(J/V 4> ^ D^Df*) + a(J/V </> ^ i?,-*^^^), 



(Dll) 

(D12) 

(D13) 
(D14) 



aK = 4x a{J/ip K DsD) + a{J/i) K D*D) + a{J/i) K D^D*) + a{J/i) K D*D*) , (D15) 



A 






£)0{*) £)0{*) 


D°D°* 


£)+(*) £)-(*) 


D+D-* 


Dt'*^D7'*^ 


DfD-* 




DsD* 


DtD 




1 


2 































1 


2 





























n" 








1/2 


1 


1/2 


1 

























1/2 


1 


1/2 


1 

















V 








1/2 


1 


1/2 


1 

















4> 




















1 


2 



















n 





n 











1 


1 


1 


K" 


























i 


i 


i 




























1 


1 


1 



TABLE IV: Multiplicity coefficients for tlie processes J/tp A DD. Notice tliat D^*'> D^*'> = DD, D*D*. 
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